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We study the stability of the fixed-point solution of an array of mutually coupled logistic maps, 
focusing on the influence of the delay times, Ty, of the interaction between the ith and jth maps. 
Two of us recently reported [Phys. Rev. Lett. 94, 134102 (2005)] that if Ty are random enough 
the array synchronizes in a spatially homogeneous steady state. Here we study this behavior by 
comparing the dynamics of a map of an array of N delayed-coupled maps with the dynamics of 
a map with TV self-feedback delayed loops. If N is sufficiently large, the dynamics of a map of 
the array is similar to the dynamics of a map with self-feedback loops with the same delay times. 
Several delayed loops stabilize the fixed point, when the delays are not the same; however, the 
distribution of delays plays a key role: if the delays are all odd a periodic orbit (and not the fixed 
point) is stabilized. We present a linear stability analysis and apply some mathematical theorems 
that explain the numerical results. 
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I. INTRODUCTION 

Cooperative behavior arises in many fields of science 
and classical examples include the onset of rhythmic ac- 
tivity in the brain, the flashing on and off in unison of 
populations of fireflies, and the emission of chirps by 
populations of crickets Of important practical ap- 
plications are the synchronization of laser arrays and 
Josephson junctions. Coupled map lattices (CMLs) |2( 
are excellent tools for understanding the mechanisms of 
emergency of synchrony in complex systems composed of 
interacting nonlinear units, because, by simplifying the 
dynamics of the individual units, CMLs allow the simu- 
lation of large ensembles of coupled units. 

The effects of time delays arising from the finite prop- 
agation time of signals have received considerable atten- 
tion. Classical examples are the Mackey-Glass model in 
physiology Q (that describes anomalies in the regenera- 
tion of white blood cells due to the finite time of propaga- 
tion of chemical substances in the blood), and the Ikeda 
model Q in optics (that accounts for the finite velocity 
of light in optical bistable devices). 

Three common consequences of time-delays are mul- 
tistability, which typically arises for delays longer than 
the intrinsic oscillation period [5|,|(j, chaotic dynamics, 
which arises for strong coupling and/or long delays, and 
oscillation death, which refers to the existence of sta- 
bility islands in the parameter space (coupling strength, 
delay time) where the amplitude of coupled limit-cycle 
oscillators is zero 0. It is also well-known that time- 
delayed feedback can stabilize unstable orbits embedded 
in chaotic attractors Q and enhance the coherence of 
chaotic 0, 0| and stochastic motions 0] . 

Most studies of delayed coupling have considered uni- 
form delays, i.e., the interactions between the different 
elements of a network occur all with the same delay time 



(instantaneous coupling is a particular case of "fixed de- 
lay coupling"). To the best of our knowledge, the first 
study of a system of mutually coupled units interacting 
with different, randomly chosen delay times was don e by 
Otsuka and Chern in the early 90 's. In Ref. [ij 
an array of semiconductor lasers with incoherent optical 
coupling was studied numerically, and it was shown that 
different dynamic regimes occur, including synchroniza- 
tion, clustering and steady-state behavior, depending on 
the average delay and the delay distribution. Distributed 
(or random) delays have been the object of recent atten- 
tion since several authors have reported that non-uniform 
delays can have a stabilizing effect. Atay et al. [13} stud- 
ied ensembles of limit-cycle oscillators and showed that 
distributed delays can enlarge the stability islands where 
oscillator death occurs. Huber and Tsimring |14| studied 
networks of globally coupled, noise-activated, bistable el- 
ements and found that increasing the non-uniformity of 
the delays enhanced the stability of the trivial equilib- 
rium. Eurich et al. [15| showed that distributed delays 
increase the stability of predator-prey systems including 
two-species systems, food chains, and food webs. 

Two of us recently studied an array of logistic maps 
coupled with randomly distributed delay times |16| , 

JV 

Xi(t + 1) = f[xi(t)] + ^ £ - T ^')] ~ /M*)]) > 

J'=l 

where t > is an integer- valued time index, i — 1, . . . , N 
is a space index, f(x) = ax(l — x) is the logistic map 
(a G (0,4]), e is the coupling strength (e 6 [0,1]) and 
Tij > is an integer that represents the delay time in 
the interaction between the ith and jth maps. For ry 
random enough the array synchronizes in the spatially 
homogeneous steady-state, Xi(t) = xo for all i, where xo 
is the nontrivial fixed point (xq = 1 — 1/a). This synchro- 
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nization behavior is in contrast with the synchronization 
with fixed and distant-dependent delays. For fixed de- 
lays (jij — To V i, j) the array synchronizes in a spa- 
tially homogeneous time-dependent state, Xi(t) = x(t) 
V i, t, where the dynamics of an element of the array 
is either periodic or chaotic depending on tq 17] . For 
distant-dependent delays (ry = k\i — j\ where k is the 
inverse of the velocity of transmission of information) a 
one-dimensional linear array synchronizes in a state in 
which the elements of the array evolve along a periodic 
orbit of the uncoupled map (i.e., Xi(t) is a solution of 
Xi(t + 1) = f[xi(t)]), while the spatial correlation along 
the array is such that Xi(t) — Xj(t — Tij) V i, j (i.e., a map 
sees all other maps in his present, current, state) [Isj . 

In Ref.^| the stabilization of the fixed-point solution 
due to random interaction delay times was interpreted as 
a "discrete time" version of the control method for sta- 
bilizing a fixed point recently proposed by Ahlborn and 
Parlitz [lflj . In Ref.^j| the fixed point of a dynamical 
system x — f(x) was stabilized with the addition of sev- 
eral feedback terms that satisfy: (i) the feedback terms 
vanish in the steady state and (ii) the delay times are not 
an integer multiple of each other (with these conditions 
the control terms vanish only at the fixed points and not 
at the periodic orbits). Simulations of a single logistic 
map with N self-feedback time-delayed loops, 



x(t + 1) = f[x(t)} + (f[x(t - r 3 )] - f[x(t)]) , (2) 



N 

show that several terms with different delays lead to the 
stabilization of the fixed point after transients. 

The behavior of a single unit often helps understand- 
ing the behavior of an ensemble of coupled units, and in 
particular the "chaos suppression by random delays" in 
an ensemble of coupled logistic maps can be interpreted 
in terms of the suppression of chaos and the stabilization 
of the fixed point in a single logistic map with several 
delayed self-feedback loops. The aim of this paper is to 
further investigate this point, by comparing the dynam- 
ics of an element Xi of an array of globally coupled N 
logistic maps [Eq. QJ] with the dynamics of a logistic 
map with N self- feedback loops [Eq. J5J]. 

This paper is organized as follows. Section II presents 
a linear stability analysis of the fixed point solution of 
Eq. (J2J and discusses the stability in the parameter space 
(local nonlinearity, a, feedback strength, e). We find an 
analytical (sufficient) instability condition, Eq. l|15[l . that 
holds for large a and low e, regardless of the number of 
feedback terms and/or the values of the delay times. We 
find a second (sufficient) instability condition, Ea. (|21|l . 
that holds for one feedback loop, independently of the de- 
lay time. Section III presents a comparison of the dynam- 
ics of a logistic map with TV self- feedback delayed loops, 
with the dynamics of a map of an array of N delayed- 
coupled logistic maps. The numerical simulations show 
that if N is sufficiently large, the dynamics of a map of 
the array is remarkably similar to the dynamics of a map 



with N feedback loops and the same delay times. The 
similarities are explored by analyzing the regions in the 
parameter space (a, e) where the fixed-point solution is 
stable, and comparing with the analytic results of Sec. II. 
We also present bifurcation diagrams that demonstrate 
similar types of instability scenarios. Section IV presents 
an interpretation of these results based on the analogy 
between globally coupled maps and a single map with 
external driving, studied by Parravano and Cosenza in 
Refs.pol l2*l) . Section V presents a summary and the 
conclusions. 



II. STABILITY ANALYSIS 

To analyze the stability of the nontrivial fixed point 
solution of Eq. (2), x — f(xo) (for the logistic map xq = 
1 — 1/a), we define a new set of variables, 



y m (t) =x(t- m), 



(3) 



with m — 0, ...,M and M = max(rj), that describe the 
present and past state of the map. We can re- write Eq.J2J) 
in terms of these new variables as: 



y m -i(t) ifrn^O, 
(l-e)/bfo(t)] + £££i/[»r«] ifm = 0> ' 



The fixed-point solution is 

2/o (*) =x ;...;y M (t) = x . (5) 
To study the stability of this solution we linearize, 
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Sy m (t + !) = /] A mn 5y n (t), 



(6) 
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where 



A = 
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(8) 



where l n is the number of times the value r = n appears 
in the sequence ri, T2,... t^: 'Y^— 1 l n — N. The term 
cxq accounts for the instantaneous feedback loops. Notice 
that some of the a n coefficients will be zero (a n = if 
Tj =/= n Vi); however, the coefficient corresponding to the 
maximum delay, M — T maxi is different from zero and is 
given by 



a M 



lMjjf'(x ). 



(9) 
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The next step for the derivation of analytic stability con- 
ditions is the study the eigenvalues A^ (with i = . . . M) 
of the matrix A. The Gershgorin theorem [22j states that 
all eigenvalues of a complex square matrix are located in 
a set of disks centered at the diagonal elements an with 
radius equal to the sum of the norms of the other ele- 
ments on the same row: 



|Aj -a u \ < ^ \a>ij\ 



(10) 



For i ^ Eq. QlOfl gives \Xi\ < 1 and for i = gives 

\X i -(l- e)f'(x )-a \ < e|/'(aro)| - |«o|, (H) 



where we used Eq. (JSJ, J2j=o h = ^ anc ^ J2j=o a j ~ 
ef'(xo). Therefore, the eigenvalues are in the region of 
the complex plane defined by the two disks: 



|A| < 1 

|A-(l-e)/'(x )-a | < e|/'(x )| 



(12) 
(13) 



From here we derive a sufficient stability condition and 
a sufficient instability condition. If the disc of radius 
e | /' (xo ) | — | ceo I centered at (1 — e)f (xq ) + a is completely 
inside the disc of radius 1 centered at 0, then all the eigen- 
values will have |A| < 1. Therefore, a sufficient stability 
condition is \ (1 — e) f (x ) + a \ + e\f (xq) \ — \ao\ < 1, and 
taking into account Eq. JSJ all the a n and f'(xo) have 
the same sign, the sufficient stability condition read as 



\f(x )\ < 1. 



(14) 



This stability condition is trivial because if Eq. l|14fl holds, 
then the fixed point of the "solitary map" (the map with- 
out feedback loops, e = 0) is stable. 

Let us consider the region where \f (x )\ > 1 (for the 
logistic map \f (x )\ > 1 for a > 3). In this region we 
have the following sufficient instability condition: if 

|(l-e)/ / (x )|-e|/'(a ;o )| + 2|a | > 1, (15) 

then the disc centered at (1 — e)f'(xo) + a of radius 
e|/'(xo)| — | ao I is completely outside of the disc centered 
at zero of radius 1, and therefore there is at least one 
eigenvalue with |A| > 1. We remark that this instabil- 
ity condition holds, regardless of the number of feedback 
terms, and regardless the values of the delay times. 

The stability and instability regions are displayed in 
Fig. 1(a). The trivial stability condition holds for a < 3, 
the instability condition holds for large a and small e 
(there is a second instability region, in the corner of large 
a and e which is discussed below). For comparison, we 
show in Figs. l(b)-l(d) the stability regions calculated 
from numerical simulations of Eq.Q with one delayed 
feedback loop (and different delay times), which agree 
well with the analytic predictions. 





(d) 



FIG. 1: (a) Analytically calculated stability regions for the 
logistic map with one delay term. The trivial stability region 
defined by Eq. 11411 is displayed in black (blue online), the 
instability regions defined by Eqs. (11511 and 12 U are displayed 
in dark grey (red online), (b)-(d) Stability regions calculated 
numerically, by simulation of Eq. JSJ with one delay term and 
different initial conditions. The parameter regions where the 
fixed point is stable for all initial conditions are displayed in 
black (blue online). The black solid lines indicate the borders 
of the instability regions defined by Eqs. 11511 and 12 U . Near 
these boundaries there is sensitivity to the initial conditions: 
some trajectories evolve towards the fixed point, while others 
evolve to periodic or chaotic orbits. The delay times are (b) 
ri = 1, (c) ri = 2, (d) n = 3. 



Additional information can be obtained by calculating 
explicitly the eigenvalues of A. The roots of the charac- 
teristic equation 



/(l - e)f'(x ) + a - A cvi a 2 ... a M -i a M \ 



dot 
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1 -A 








V 



0. 



... -A 

... 1 -A/ 

(16) 

can be written in terms of the determinant of two M x M 
matrices: 



[(1 - e)f'(xo) + a - A] dct C - dct B = 0, (17) 



where 



C = 



/—A 
' 1 -A 




V 
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-A 
1 -A/ 



(18) 
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and 



B = 



fOL\ CC2 

1 -A 




V 



OiM-l Ct M \ 





-A 
1 -A/ 



(19) 



We calculated the determinant of each matrix recursively 
(details are presented in the appendix) obtaining: 

M 

A M+1 - [(1 - e)f(xo) + ao]A M -J2^X M ~ j = 0.(20) 

The roots of Eq.® satisfy ]jfi Q = \a M \. Thus, if 

\a M \ = \l M ef'(x )/N\ > 1 (21) 

at least one eigenvalue has |A| > 1, i.e., this gives another 
analytic (sufficient) instability condition. For Im/N = 1 
Ea. (|21|l holds in the right-bottom corner of Fig. 1(a) 
(large a , large e); for Im/N < 1 Ea. H21fl is not satisfied 
in the parameter region of interest (a € [0,4], e € [0,1]). 
Notice that if Im/N — 1, then the map has only one 
feedback loop (because the multiplicity of the feedback 
terms with maximum delay is equal to the total num- 
ber of feedback terms); therefore, Ea. (|21fl indicates that 
when a logistic map has a single feedback term, the fixed- 
point solution can not be stable in the parameter region 
(large a, large e), regardless of the delay time. In other 
words, a single feedback term can not stabilize the fixed 
point in the (large a, large e) parameter region. Numeri- 
cal simulations of Eq. J5J with N = 1 and different delays 
confirm these analytical predictions, see Figs. l(b)-(d). 

Further analytical insight can be gained by considering 
two special cases: all-even and all-odds delays. First, let 
us show that if the delays are all even (and therefore, M 
is even), A = —1 is a solution of Eq. lj2U|) when f'(xo) = 
— 1, i.e., at the border of the stability region, Ea. (ll4|l . 
Substituting A = — 1 in Ea. (|20|) and taking into account 
that M — j is even (since in the sum only terms with j 
even arc different from zero) gives 



1 - [(1 - e)f'(xo) + ao] 



M 

£ 

3=1 



(22) 



Using J2jLo a j ~ e f'( x o) we obtain f'(xo) = — 1. There- 
fore, when the delays are all-even there is an eigenvalue 
A = — 1 if and only if a = 3, regardless of e. 

Next, let's see what happens if the delay times are all 
odd, therefore, M is odd and, in addition, ao = 0. Taking 
into account that M — j is even (since in the sum only 
terms with j odd are different from zero) for A = — 1 Eq. 
(|2"U|l gives 
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FIG. 2: Eigenvalues in the complex plane for a map with 
N = 3 feedback loops with all-odd delays (left column: tj =1, 
3, 5) and all-even delays (right column: Tj —2, 4, 6). e = 0.5 
and (a), (e) a = 2.5; (b), (f) a = 2.7; (c), (g) a = 3.0 and (d), 
(h) a = 3.8. The circles indicate the Gershgorin disks. 



Using Y,j=i a j = e/'(a;o) we obtain f(x ) = -l/(l-2e). 
For the logistic map this gives a — (3 — 4e)/(l — 2e) a 
condition which is satisfied for e 6 [0, 1] only for values 
of a inside the stability region a < 3. 

The above analysis allows as to draw some additional 
conclusions about the stability of the fixed point in the 
special cases of all-even and all-odd delays. For all-even 
delays an eigenvalue is real and negative and equal to — 1 
for a = 3. For larger a this eigenvalue can in principle 
become A < — 1 rendering the fixed point unstable due 
to a period-doubling bifurcation. For all-odd delays this 
instability scenario is not possible, as A = —1 in a pa- 
rameter region where we know that all eigenvalues must 
have |A| < 1. 

We verified these predictions by calculating numeri- 
cally the eigenvalues of A. Figure 2 displays results for 
all-odd and all-even delays, varying a while keeping e 
constant. It can be observed that for all-even delays one 
real eigenvalue becomes less than —1 for a > 3; for all- 
odd delays a pair of complex-conjugate eigenvalues have 
modulus greater than 1 for a > 3.8. 
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III. NUMERICAL RESULTS 

In this section we compare the dynamics of a logis- 
tic map with N self- feedback loops, Eq. with the 
dynamics of a map of an array of N globally delayed- 
coupled logistic maps, Eq. (J). 



We consider Gaussian distributed delays: Tj. 



To 



Near(c£), where c is a parameter that allows varying the 
width of the delay distribution (for c = the delays are 
all equal, = To, for c =/= the delays are distributed 
around To); £ is Gaussian distributed with zero mean 
and standard deviation one; Near denotes the nearest 
integer (we use Near instead of Int to have a distribution 
that is symmetric with respect to To; however, the results 
are largely independent of the precise form of the delay 
distribution). Depending on to and c the distribution of 
delays has to be truncated to avoid negative delays. 

We begin by showing that if TV is large enough, the 
parameter region where the fixed point is stable for the 
array is remarkably similar to the parameter region where 
the fixed point is stable for the map with feedback loops. 
Figure 3 displays results for three values of N, the up- 
per row shows the stability region of the homogeneous 
steady-state solution of the array [xi(t) = Xj(t) — xq V i, 
j, t], while the lower row shows the stability region of the 
fixed point solution the map with feedback loops. The 
delays of the self-feedback terms, Tj with j — 1 . . . N, 
were taken equal to the delay times of the interaction of 
the ith and jth maps of the array; this gives N sets of 



delay times (t, 



Tij with i 



1 . . . TV). In the lower row 



of Fig. 3, the parameter region where the fixed point 
is stable for all sets of delays is displayed in black (blue 
online), and the regions where is unstable for all sets of 
delays are displayed in dark gray (red online). Outside 
the trivial stability region (a < 3) and outside the insta- 
bility region defined by the sufficient condition Eq. 1)150. 
if TV is small the fixed point of the map with self- feedback 
loops can be stable or unstable depending on Tj (i.e., the 
fixed point can be stable for the ith set of delays and not 
for the feth set of delays); however, if TV is sufficiently 
large the stability of the fixed point is the same for all 
sets of delays (there is sensitivity to the precise values 
of Tj near the boundaries of the instability regions). It 
can be observed that for both, the single map and the 
array, the fixed point is unstable in the left-bottom cor- 
ner (large a, low e), in agreement with the results of the 
previous section (where we showed that in this region the 
sufficient instability condition, Eq. l|15|l . holds). 

The stability of the fixed point depends on the distri- 
bution of delays, and again, the similarities between a 
map of the array and a map with self-feedback loops, for 
TV large enough, are remarkable. As the width of the de- 
lay distribution, c, increases, the parameter region where 
the fixed point is stable grows (see Fig. 4), and this oc- 
curs for both, the array and the map with self-feedback 
loops. 

For the array of coupled maps, the parameter that 
quantifies the influence of the delays is not the mean 




(c) 
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FIG. 3: Stability region of the fixed point solution of an array 
of TV maps (upper row) and of a map with TV feedback loops 
(lower row). The delays are Gaussian distributed with to = 3, 
c = 1. The dynamics of the map with feedback loops was 
simulated for various sets of delays Tj: the fixed point was 
found to be stable for all sets in the black region (blue online) 
and unstable for all sets in the dark gray regions (red online). 
The solid line indicates the borders of the instability region 
defined by Eq. (GJ. (a) TV — 10, (b) TV = 20, (c) TV = 100. 
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FIG. 4: Influence of the width of the delay distribution on 
the stability of the fixed point solution of an array of TV maps 
(upper row) and of a map with TV feedback loops (middle 
row). The delays (lower row) are Gaussian distributed with 
to = 3, (a) c = 0.5, (b) c = 1 , (c) c = 2. TV = 100. 



delay, < t >, or the standard deviation of the distri- 
bution, D T , but is the normalized disorder parameter, 
c* = D T / < t > [3 (< t >= T , D T = c, and c* = c/t 
if the Gaussian distribution is not truncated). Figure 5 
shows that this is also the case for the map with self- 
feedback loops, as it can be noticed that the stability re- 
gion of the fixed point is the same for distributions that 
have different < r > and D T , but the same normalized 
width, c*. 
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FIG. 5: Influence of the normalized disorder parameter on 
the stability of the fixed point solution of an array of N maps 
(upper row) and of a map with N feedback loops (middle 
row). The delays (lower row) are Gaussian distributed with 
different mean and different standard deviation but with the 
same normalized width, c* ~ c/ro (a) To = 3, c = 1, (b) 
T = 6, c = 2 , (c) r = 10, c = 3.33. N = 100. 



Figure 6(a) displays the stability region of the fixed 
point solution for the array (upper row) and for the map 
with self-feedback loops (lower row), in the parameter 
space [to (~< r >), c (~ D T )\. It can be observed that 
the fixed point is stable if the delays are random enough 
(i.e., if the width of the distribution, c, is larger than a 
certain value that increases with To). When the stability 
region is plotted vs. the normalized width, c* [Fig. 6(b)], 
it can be observed than the value of c* above which the 
fixed point is stable is independent of r [but depends on 
e, as shown in Fig. 6(c)]. This occurs for both, the map 
of the array and the map with self- feedback loops. 

While the above observation of enhanced fixed-point 
stability with increasing randomness of delays is generic, 
independent of the precise form of the delay distribution, 
there is an exception which is the case of all-even delays. 
For all-even delays the fixed point of the map with self- 
feedback terms is stable only in the region defined by the 
sufficient trivial stability condition, Ea.ip.4j). Because the 
fixed point becomes unstable due to a period-doubling 
bifurcation when one real eigenvalue becomes A < — 1 (as 
discussed in the previous section) , all-even delays tend to 
stabilize an orbit of period 2. The same effect is observed 
in the array of N logistic maps: for all even delays the 
fixed point is stable only in the trivial region a < 3. 

If N is sufficiently large, a map with N self-feedback 
loops and a map of an array of N maps follow very similar 
instabilities scenarios when e or a are varied (see below for 
a discussion of the limit in which the similarities are not 
only qualitative but also quantitative). As an example, 
Figs. 7-9 display bifurcation diagrams for varying e while 
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FIG. 6: (a) Stability region of the fixed point solution of an 
array of N maps (upper row) and of a map with N feed- 
back loops (lower row) in the parameter space (mean delay, 
standard deviation), (b) Same as (a) but plotted vs. the nor- 
malized disorder parameter, see text, (c) Same as (a) but in 
parameter space (e, disorder parameter). N = 100, a — 4, in 
(a),(b) e = 1; in (c) To = 5. 



keeping a fixed. The delays are "mixed" (even and odd) 
in Fig. 7, all-odd in Fig. 8, and all-even in Fig. 9. The 
parameters correspond to a scan of e along the horizontal 
axis of Fig. 4(b): for "mixed" delays the fixed point is 
stable in a range of e for both, the map with feedback 
loops and the array. Figures 7(a)-9(a) [7(b)-9(b)] display 
the time evolution of one element of the array, i = 1 (i = 
2), by plotting 100 consecutive interactions of X\ [x<x] 
after transients die away vs. e. Figures 7(c)-9(c) display 
the array configuration at time t — T (large enough to let 
transients die away) . Figures 7(d)-9(d) [7(e)-9(e)] display 
the time evolution of the map with feedback loops with 
delays Tj — T±j (tj — T2j), by plotting 100 consecutive 
interactions of x (after transients die away) vs. e. 

Above a certain coupling strength the array synchro- 
nizes in a single cluster that displays either steady-state 
or time-dependent dynamics [in Figs. 7(c)-9(c) there is a 
single cloud of points for e >« 0.5]; below this coupling- 
strength the array splits into clusters and the elements 
of each cluster evolve along similar time-dependent or- 
bits (the bifurcation diagrams for the elements i = 1 and 
i — 2 of the array are similar, even for low e). 

For a map with self-feedback loops with "mixed" de- 
lays, the fixed point is stabilized for increasing e after a 
period-doubling bifurcation [Fig. 7(d), (e)]. In contrast, 
for feedback loops with all-odd delays the fixed point is 
stabilized after a Hopf bifurcation [Fig. 8(d), (e)]. For 
all-even delays the fixed point is not stable for any e [but 
the period-two orbit is stable in a certain range of e, 
Fig. 9(d), (e)]. These results are in agreement with the 
analysis of the previous section, where we found that for 
all-odd delays the fixed point changes stability when a 
pair of complex eigenvalues cross the unit circle, and for 
all-even delays the fixed point changes stability when a 
real eigenvalue crosses the unit circle at A = — 1. The 
bifurcation diagrams of the i = 1 and i = 2 maps of the 
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FIG. 7: Bifurcation diagrams for increasing e. (a) i = 1 
map of the array, (b) i — 2 map of the array, (c) Array 
configuration, Xi with i = 1, . . . , N at time t = T. (d) A map 
with feedback loops with delays Tj = Tij. (e) A map with 
feedback loops with delays Tj = Tiy. a = 4, N = 100, the 
delays are Gaussian distributed with to = 3, c = 1. 
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FIG. 8: As Fig. 7 but with all-odd delays. 



array display similar features [Figs. 7-9(a), 7-9(b)]: the 
fixed point is stable in a range of e for "mixed" and all- 
odd delays, while the period-two orbit is stable in a range 
of e in the case of all-even delays. 

Exponentially distributed delays [r^ = To + Int(c£), 
where £ is exponentially distributed, positive, with unit 
mean] yield similar bifurcation diagrams, shown in Fig. 
10. Furthermore, the instability scenario for fixed e and 
increasing a (i.e., a scan along a vertical line in Fig. 4) 
is also very similar, as shown in Fig. 11. 
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FIG. 10: As Fig. 7 but with exponentially distributed delays 
(ro = l, c=l). 



IV. DISCUSSION 

The similarities between a map of the array and a map 
with self-feedback loops can be interpreted in the frame- 
work of the analogy between globally coupled maps (with 
instantaneous coupling), and a single map subjected to 
an external drive, studied by Cosenza and Parravano in 
Rcfs. 20, 211. The authors considered 



Xi {t + 1) = (1 - e)f[xi{t)} + eH( Xl (t) . . . x N (t)), (24) 

where H is a global coupling function that is invariant to 
argument permutations [H(x\, . . . , Xi, . . . , Xj, . . . xn) — 
H(xi, . . . , Xj, . . . , Xi, . . . , xn) V i and j], and showed that 
the clustering behavior of the array can be analyzed 
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FIG. 11: As Fig. 7 but varying a while keeping e = 0.3 fixed. 
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£ 

FIG. 12: (a) Mean field for the i — 1 map of the array vs. e. 
(b) Mean field for each map of the array, Hi with i = 1 . . . N, 
vs. e. (c) Driving term of the map with feedback loops, F, 
vs. e. N = 100, a = 4, the delays are the same as in Fig. 7. 



through the analogy with the driven map, 

x(t+l) = (l-e)f[x(t)] + eF(t), (25) 

where F(t) is a external forcing (assumed to be periodic). 
The analogy holds because in Ea. (|24l) all the elements 
of the array are affected by the coupling function H in 
exactly the same way at all times, and therefore the be- 
havior of any element of the array is equivalent to the 
behavior of the driven map, Eq. i|25|) . 

To analyze whether this analogy can be extended to 
the case of delayed-coupling, we calculated the mean field 
coupling term at site i of the array, 

i N 

ffiW = ^EW-4 (26) 
i=i 

and compared with the driving term of the map, 
1 N 

F(i) = ^Y,f\- x{t - T ^- ( 27 ) 

3=1 

We found that for N large, Hi(t) is nearly the same for 
all the elements of the array, regardless of the array dy- 
namics, and its time evolution is similar to that of the 
driving term of the map, F. As an example, Fig. 12 dis- 
plays the mean field coupling term and the driving term, 
for parameters corresponding to the bifurcation diagrams 
shown in Fig. 7. In Fig. 12(a) we plot 100 consecutive 
values of the mean field at one element of the array, H\ 
(after transients die away). Figure 12(b) displays the 
mean field at all elements, Hi with i = 1 . . . N, at time 
t = T (large enough to let transients die away). It can 
be seen that Hi ~ Hj even for low values of e. Figure 
12(c) displays 100 consecutive values of the driving term 
of the map with feedback loops, F(t) (after transients die 



away), and it is observed that above a certain value of e 
(e w 0.2) F ~ H v 

We speculate that the analogy with the single map has 
its roots in the fact that the elements of the array display 
similar temporal variation, i.e., they evolve along equal 
(or similar) orbits, even when the array splits into clus- 
ters. Therefore, a map of the array "perceives" signals 
coming from other maps as nearly indistinguishable from 
signals coming from self-feedback loops. This can also be 
thought as an ergodic property of the dynamics, since 
the average over the ensemble, Hi, is nearly equal to the 
average over time, F. 

The analogy holds even if the delays are all equal, 
c = 0, as shown in Fig. 13. Above a certain cou- 
pling strength (e « 0.45) the array synchronizes in-phase, 
Xi(t) = xj(t) V i, j, and therefore the analogy is math- 
ematically trivial since in the synchronization manifold 
the evolution equation of one element of the array and the 
evolution equation for a map with a single self-feedback 
loop are exactly the same. However, multistability in 
the delayed map (i.e., the coexistence of different stable 
orbits) might lead to competition phenomena in the ar- 
ray of coupled maps, as different elements might tend 
to evolve along different orbits, depending on the initial 
conditions. The investigation of this type of dynamics is 
the object of future work. 



V. SUMMARY AND CONCLUSIONS 

We investigated the dynamics of mutually coupled lo- 
gistic maps focusing on the influence of the delay times of 
the interactions between the maps, and comparing with 
the dynamics of a map with several time-delayed self- 
feedback loops. By using some mathematical tools such 
as the Gershgorin theorem we derived analytic stability 
and instability conditions for the fixed point solution of 
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FIG. 13: As Fig. 7 but with fixed delays (r = 1, c = 0). 



the map with feedback loops. We found that the stabi- 
lization of the array in the fixed-point solution can be 
well understood in terms of the dynamics of the map. 
Specifically, for randomly distributed delay times, if N 
and e are large enough the fixed-point is stable for both, 
the map with self-feedback loops and the array. Also, 
if the delay times are all even, for both, the single map 
with N delayed loops and the array of N delayed-coupled 
maps, we observed that the stability region of the fixed- 
point is reduced to the "trivial" region (a < 3) regardless 
of the coupling strength. The results presented here pro- 
vide another example of an ensemble of mutually coupled 
interacting units, where understanding the dynamics of 
a single unit with self-feedback loops is relevant for un- 
derstanding the macroscopic behavior of the ensemble. 
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VII. APPENDIX 



In this appendix we demonstrate Eq. I|2(J|) . We start 
from Eq. ljT7|) . 



[(1 - e)f(x ) + a - A] det C - det B = 0, (28) 
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(29) 



The second determinant is zero because the elements of 
the first row are all zero. We obtain 



detC(M) = -A det C(M - 1) 



(-A) 



M 



(30) 



The determinant of matrix B can also be calculated re- 
cursively: 



det B = det 



1 -A 
1 -A 





V o o o 



a.\ det 
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Qfftf-l a M\ 




-A 
1 -A/ 



(31) 



It can be noticed that the first matrix is C(M — 1) while 
the determinant of the second matrix can be calculated 
recursively, 

det B(M) = ai det C(M - 1) - det B(M - 1), (32) 
where 

/ a 2 ••• QfM-i o>m\ 



B(M — 1) = 




V 



-A 

1 -a; 



(33) 
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Substituting in (|28(l that can be simplified to 

[(1 -e)/'(x ) + a - A] dct C(M) — [a\ det C(M - 1) 
-a 2 det C(M - 2) + det B(M - 2)] = (34) 



Using Eqs.lOJJ and(j23 we obtain 



[(1 - e)f(x ) + a - A](-A) M -« 1 (-A) M - 1 

+a 2 (-\) M - 2 - ••■-0 (35) 

which gives 

(-X) M+1 + [(l-e)f(x ) + a ](-X) M 

M 

+ J2(-l) j a j (-X) M - j =0, (36) 
i=i 



M 



A M+1 + [(1 - e)/'(x ) + a ]A M + ^ „,,V U ' = 0(37) 
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